#include <complex.h>
#include <stdio.h>
#include <stdlib.h>
#include <math.h>

#include "constants.h" 

#include "BiCGSTAB_precond_CRS.h"
#include "QMR_precond_CRS.h"



void MaxwellVFEM(double complex *A_EM, double complex *phi, 
		 double complex *Jc, double complex *rho, 
		 double complex omega, double complex *epsilon, double *mu, 
		 int N_node, int N_edge, 
		 int N_triangle, int N_triangle_interface, int N_element, 
		 double *Node, int *Edge, int *Triangle, int *Element, 
		 double *Length_of_Edge,  int *Triangle_Edge, 
		 int *Element_Edge, int *Triangle_Element);

/* void MaxwellVFEM(double complex *A_EM, double complex *phi, double complex *Jc, double complex *rho,  */
/* 			   double complex omega, double complex *epsilon, double *mu,  */
/* 			   int N_node, int N_edge, int N_element,  */
/* 			   int N_point, int N_line, int N_triangle,  */
/* 			   int N_line_interface, int N_triangle_interface,  */
/* 			   double *Node, int *Edge, int *Element,  */
/* 			   int *Element_Edge, int *Triangle_Edge, int *Line_Edge,  */
/* 			   int *Point, int *Line, int *Triangle,  */
/* 			   int *Triangle_Element,  */
/* 			   double *Length_of_Edge); */


